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Abstract 

The asymptotic beliavior of three-body scattering wave functions in config- 
uration space is studied by considering a model equation that has the same 
asymptotic form as the Faddeev equations. Boundary conditions for the wave 
function are derived, and their vahdity is verified by numerical calculations. 
It is shown that these boundary conditions for the partial differential equation 
can be used to obtain accurate numerical solutions for the wave function. 



I. INTRODUCTION 

In a previous paper [|T| (hereafter referred to as I), we studied the asymptotic form of 
a three-body wave function, which results in the propagation of three free particles from 
various types of sources. The goal was to determine the boundary conditions appropriate 
for the three-body scattering equations in configuration space. In particular, we want to 
establish the values of 
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for which the leading asymptotic form of an outgoing wave would be valid and to investi- 
gate the correction terms to that form. For a complete discussion of three-body scattering 
in configuration space see Ref. 0, which also contains references to earlier work on this 
problem. 

In I two sources were studied. The first was a localized source corresponding to the 
elastic-scattering driving term in the three-body Faddeev equation and is determined by 
the overlap of a two-body force in one pair and a two-body bound-state wave function in 
another pair. The other was an extended model source that mimics the real source term in 
the Faddeev equation including the breakup process in the Faddeev amplitude. The latter 
source reaches far out in the distance y between one particle and the center of mass of 
the other two particles; it decreases only as 0{y^^^'^). However, the presence of the pair 
interaction limits the extent in the distance x between the other two particles. As expected 
in the case of the extended source, the leading form is reached only at a much larger radius 
than the localized source, specifically when x is small and y is large. 

By inverting the free propagator, one can determine the propagating wave function using 
a partial differential equation with the given source terms. We established suitable boundary 
conditions that could be used to solve this problem efficiently. A matching radius of about 
100 fm was found to be sufficient. 

In this article we extend our previous study of the extended source to allow one pair 
interaction to be present while the three particles propagate from the given sources, which 
is exactly what happens in the Faddeev formulation. We use the same notation as in I. In 
Sec. II we evaluate the three-body Green's function including one pair interaction, apply 
it to the extended source, and study the asymptotic behavior in the two- and three-body 
fragmentation channels. In Sec. Ill we solve the related partial differential equation as an 
exercise for applying that technique (in a forthcoming article) to the Faddeev equation itself. 
Finally, we conclude in Sec. IV. 
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II. THREE-BODY PROPAGATION FROM GIVEN SOURCES 
WITH A PAIR INTERACTION 



In order to avoid unnecessary complications, we restricted our considerations in I to three 
identical bosons interacting by spin-independent s-wave pairwise interactions in a state with 
total angular momentum of zero. For this case the Faddeev equation for the channel with 
incident wave 

= ^ + G{E)VP^ , (1) 
where P = + P~ is the permutation operator, reads in explicit notation as 

/■OO /"OO 

i^{x,y) = (f){x,y) + dx' dy' g{x,y; x',y') Q{x' ,y') . (2) 
JO Jo 

We have introduced in Eq. the reduced Faddeev amplitude ipi^x, y) and the corresponding 
reduced Green's function g{x,y;x',y'). The coordinates x and y are the standard Jacobi 
variables 

X = xi = r2 - rs (3) 
y = yi = ri - |(r2 + ra) (4) 

expressed in terms of the individual position vectors. 
The source term in Eq. (0) is given by 

/I jry 
dfi ip{x2,y2) , (5) 

where X2, 2/2, and fi in Eq. result from (|]) and (H) by the cyclical permutations, and are 
explicitly given by 



X2 = J^x^ + y'^ + xyn (6) 



y2 = ^Jf^x^ + \y^ -jxyjj , (7) 

where fi is the cosine of the angle between x and y. The total wave function is given by 
the sum of the three Faddeev amplitudes. In Eq. (|l]) is one of the Faddeev amplitudes, 



the other two are generated by cychcal permutations of the particles, P"'"\l/ and P~\I/, and 
they appear in the source term as P"^. The pair interaction is V{x) and this interaction 
also occurs in the Green's function 

G{E) 



E-{Ho + V)+ie' 

We refer to Ref. for the general background and details on the notation. 

The source term has a short-range component arising from the elastic-scattering piece of 
the Faddeev amplitude, and a long-range component from the breakup piece of the Faddeev 
amplitude. In I we studied the effects of both components; however, for the long-range com- 
ponent we used a model source term. Using the asymptotic form of the Faddeev amplitude 
derived by the stationary phase approximation, one finds that the asymptotic form of 
the source term for three equal mass particles with total energy E is 



Q{x,y)^CVix) ^3^^ , (8) 

where fcg = mE/h^ and the constant C is given by the magnitude of the wave function in 
the asymptotic region. Therefore, to study the effects of this long-range behavior, we used 
the model source term 



xy 6 

<5Modei(a;, y) = V{x) - — ■ , (9) 

[y + yof 

with I/O = 2fm. This source term has the same asymptotic form as (|^), and we have set C 
to be unity for convenience. 

In I we neglected the final-state interaction between one pair in the propagator. As 
required by the Faddeev scheme, this will now be included. We will study the second term 
in Eq. @) with the source term replaced by the model source. Thus, we write 

poo rco 

F{x,y)= dx' dy'g{x,y;x\y')QModeiix',y') . (10) 
Jo Jo 

To simplify the numerical calculations we follow the procedure used in I, and use the 
Bargmann two-body potential 



V{x) 



;i + /3e-^^)2 ' 



where 



M 



The bound-state wave function for this potential is given by 



where k is the bound-state wave number for a two-body state with the energy e = —h'^K^/m 
and 

X + 2k 
X-2k' 

For our model calculations we use the values k, — 0.2316 fm~^ and A = 0.7 fm~^. We also 

need the two-body scattering states Ui.{x) for this potential. They are given by 



) = 1 U^^''^e'''''h{k,x) - e-'^^''^e-'''''h*{k,x)} , 



with 



1 + 



h{k, x) — 



2k -i\ 



2k + iX 



13 e 



-\x 



1 + (3 e-^^ 



and 



\ 



2k + iX 



2k -iX 



\ 



k + iK 



k — IK 



In addition, we use the arbitrary but fixed laboratory energy of the incident particle of 
-E'lab — 14 MeV and fi^/m — 41.47 MeV • fm which corresponds to the case for three nucleons. 
For this case k^ — 0.41403 fm~^. Henceforth, we set h — 1. 

For three equal-mass particles with mass m and total energy E, the three-body Green's 
function g[x, y; x', y') for the case with one pair interaction has the well-known form 



/ /N / N / 4m ■ „ sin qr.y^ \ , 



+ -/ dkuj,{x)[— ^^^(x'), 11 



where = ^Am/'i{E — e), = y4/3(A;o — and /cq = mi?. Obviously the propagation 
from the source can now proceed not only into the unbound states but also into the deuteron 
channel. Moreover, the two-body scattering states ^^(x) include the interaction V{x). It 
is the purpose of this paper to study the additional effects of V{x) in the propagator, in 
contrast to I where only the free propagator was considered. The numerical evaluation of 
the second part in (^) containing Uf.{x) requires some explanation. We were not able to 
find an analytical expression such as we found for qq in I, and had to perform the fc-integral 
numerically. Clearly at the upper end of the integral the integrand has rapid oscillations 
that make the integral difficult to evaluate numerically. Since 'U^(x) approaches sinkx as 
A; — s> oo, it appears natural to subtract the free propagator gQ{x, y; x', y') and add it back in 
a separate term. Then for A; ^ oo the integrand has a stronger fall off and, moreover, the 
path of integration can be rotated into the complex plane, similar to the procedure used in 
Ref. i. 

Let us start with the propagation in the deuteron channel: 



P / X _ 4m ( . . , e^''o^> smgoy< p 
Fd[x, y) = — — Ud[x) / dy / dx Ud[x ) QuodeAx , y 

6 Jo n„ "'0 



/„\ /"°° J ^ e^%y> singoy< 
% 

Ud{x)e''oy dy / dx Ud{x ) Quodeiix ,y 



3 Jo Jo 

4m /"°° J /Singo(2/ — 

% 

We see that this consists of a flux-conserving term 



Udix) r dy' ''""^"^^ ^'^ r dx'udix') QmoU^', y') . (12) 
Jy On -^0 



Fr{x,y) =Ud{x)e''^yfd, (13) 

with 



Am ,sm%y' p 

Id = — — / dy / dx Ud[x ) QMode\[x , y 

3 Jo q Jo 
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and a correction term 

y) ^ u,{x) r dy' ''""^"^^"^^^ H dx'u,{x') QmoU^' . y') • (14) 
6 Jy Jo 

Inserting Eq. (|^) into Eq. (|14D and performing one partial integration in the y' variable, 
one arrives easily at the following asymptotic form 

Jy/I/3kgy ^ 

- — JJ^ — Udix) - 
y^oo y^i'^ e Jo 

Clearly the long-range source behavior carries over into a corresponding long-range correction 
term in the deuteron channel. We rewrite Eq. ( [I^ ) in the form 



Fr'{x,y)^ Tj^Udix)- dx'u,{x')x'V{x'). (15) 



e 



*v ^/"^ 



Fn^, y) = u,{x) friy) , (16) 



y3/2 



and show in Fig. 1 the behavior of fd""{y) as it approaches its asymptotic value given by 

friy) ^ - - / dx'u,{x')x'V{x') . (17) 

y^oo € Jo 

To illustrate the convergence, we have normalized the plot to the asymptotic value 
/j=°"(oo) = -23.325 fm3/2. 

To illustrate the error in the elastic term that results from matching to the asymptotic 
boundary conditions at a finite distance, we plot in Fig. 2 the absolute value of fd{y) = 



Fd{x,y)/ud{x) and its asymptotic form given by Eq. ([T3|) and Eq. (15). The difference is 

less than 2% for y greater than 50 fm and less than 1% for y greater than 75 fm. 

The propagation into the unbound states u^{x) is more complicated, since 

m ^ ^ 4m 2 /■°° / . , ,e*«fc^> sing;,y< /"^ , , , ,^ ^ ^ > n /ion 
FscAx,y) = — / dkuk{x) dy / dx u^ix ) Quodeiix ,y ) (18) 

6 IT Jo Jo n. Jo 



3 1 ~~K\-^ / I 

IT Jo Jo 

can be written in the form 

4"^ 2 [''o , singly' 



FscAx,y) = — 7r~ dku^[x)e'^y dy / dx Uj^[x ) QuodeiKx ,y ) 

6 Ti Jo Jo Jo 

— / dkUf^{x) dy / dx Uf^{x ) Quodeiix ,y ) 

6 71 Jo Jy q^ Jo 

4m 2 p fy , sinh Ky' /~ , , . ,x ^ . / /^ 

— — -/ dku^{x)e y dy / dx u^,[x ) QuodeAx ,y ) 

3 7r Jko Jo K Jo 

4m 2 sinh Ky p , _^ , r ^ ' t '\ n ( ' >\ no^ 
dku,.{x) — / dye ^ / dxu^{x)QMode\{x ,y) (19) 



3 TT Jko K 

= F^^l,{x, y) + Fil(x, y) + F^l{x, y) + F^flt{x, y) . (20) 



In the third and fourth terms K = y4:/3 yk'^ — kg. 

Let us first examine the asymptotic behavior for fixed x and y approaching infinity. One 

has 



F^clt{x,y) = I ' dku,{x)e'''^yT{k) (21) 

•J 



with 



4m 2 p , singly' p , , . /^ ^ . , ,^ f^^. 

T{k) = — — - / dy / dx Uf.{x ) Quodeiix ,y). (22) 

6 n Jo n, Jo 



Am f°° J , sin Qf^y' 
% 

There is no saddle-point for this case; thus, the asymptotic form arises from the leading 
end-point contribution at A; = 0, which is easily evaluated to be 



/tt p*V of 

Fil{x,y)\,^, ~y^3'/'kfe^-/''--^u,{x)fo, (23) 
where Uq{x) = Ui.{x)/k\k=Q and Tq = T{k)/k \k=Q- We find that 



4m 2 r°° , si'^ V | ^ol/' [°° 
To = — 5--/ dy' ^= / dx'uQ{x')Quode\{x\y') . (24) 

O TT JO W i JO 



3 ^0 

This term has the same dependence on y as the correction term Eq. (^) in the deuteron 
channel. Note that Tq is given by the analytical expression Eq. ( P^j) differentiated with 
respect to k under the integral. This is obviously true for the localized source; however, 
for the extended source, one must rewrite the integral using a contour deformation before 
performing the differentiation. 

Let us now consider ^^^^^^(x, ?/) in Eq. ( [19|) for the model source term given by Eq. (P) 

(0\ 4:171 2 f^O /"OO 

Fsclt{x,y) — > — —- dku^ix) dx' u^{x')x'V{x') 

6 TT Jo Jo 



After one partial integration one finds 



(2) 2 f^o fji [°° ^^y'^/^^oV 

FscLix,y) ^ - - dkuj,ix)—J^ dx'u^{x')x'V{x')—-^ . (26) 

X fixed 



The term Fg^l^{x,y), as it stands, is less obvious in its asymptotic behavior. The contri- 
butions to the y'-integral keep growing towards the upper hmit y. Because of the factor 
exp{—Ky), only contributions from the upper end of the y' integral have to be considered. 
Again by partial integration, one easily finds 



r.(3) /X 4m 2 eV4/3fcoy roc 

^scat(a;, y) —Z, - ^ - — -TB — /. dku^ix) — 



X / dx'uk{x')x'V{x') . (27) 

<^ 

Finally, the last piece, F^^l^{x,y), can again be handled in a straightforward manner with 
the result 



(4) 4m 2 e*v^''^''o^ 
Fsclt{x,y) i -— TTo — / dku^{x) -— 



X / dx'u,^{x')x'V{x') . (28) 
Jo 

Adding equations (PB|), (P?]), and ( PBD we obtain the concise result given in Ref. 0, 

2 e«V^473fc(,y 



p(2) p(3) p(4) 



/ dku^{x)— / c/x'Mfc(x')xV(x') . (29) 

./ K J 



This can be simplified by using a technique suggested by C. Gignoux 0. Writing the 
two-body Green's function in the form 

1 2 [°° 1 

g2{x,x']z) =Ud{x) Ud{x') + - dku^{x) rr^^kW) ^ (30) 

z — e TT Jo z — k^ jm 

the integral over k occurring in Eq. (p9|) can be rewritten as 

2 /""^ m 1 

/ dkuJx) —uJx') = g2(x,x';0) +Ud(x) -Ud(x') . (31) 

Tx Jo k^ e 

Thus, we are led to the function 

roo 

u{x)= dx' g2{x,x']Q)x'V{x') , (32) 
Jo 

which obeys the inhomogeneous equation 



9 



1 rf2 



+ V{x) u{x) = —xV{x) . 



m dx"^ 

Using the explicit form for g2{x,x'] 0), one easily derives 



(33) 



u[x) — ' ~ I dx' Uq{x')x'V{x') = a 

x^oo Jo 



(34) 



where a is the scattering length defined by 6{k) approaching Koes to zero. From 

Eq. (ID and Eq. (0) follows 



u{x) = Uq{x) — X . 



(35) 



We now have the concise form 



(4) 



F^^^ -I- T?^'^^ _|_ p 
-'^scat ~r -'^scat ' -'^scat 



y-^oo y. 
X fixed 



u{x) + Ud{x) - / dx' Ud{x')x'V{x') 
e Jo 



3/2 



Altogether -Fscat(a^, has the asymptotic form 



g*V4/3 kgy 

y^oo ^3/2 
X fixed 



- u,{x)^- 3 e^-/4 To + So(a:) - x 



1 /""^ 

+'Ud(x) - / dx' Ud{x')x'V{x') 
e Jo 



(36) 



(37) 



The x-dependence is therefore built up of the zero-energy scattering state, a linear term in 
X, and the two-body bound state. The last term cancels exactly against the correction term 
Eq. ( |15|) in the deuteron channel and the total amplitude F{x,y) behaves asQ 



F{x,y) —.Ud{x)e'^'^yfd + 



y^oo 
X fixed 



7/3/2 



X] — X 



(38) 



For X outside the range of V{x), the expression in the brackets in Eq. (|38|) reduces to 

_(^_a)^3 3/Ve^"/'T^o-a. 



(39) 



^An alternate derivation of this result is given in Section 6.3 of Ref. [|[; however, the solution 
g{x) = —1 of Eq. (2.6.19) in Ref. Q is not given explicitly. 
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To verify the validity of the asymptotic term in Eq. (^) , we use Eq. ([T0|) to numerically 
evaluate F{x, y) for several values of y using only the breakup component of the Green's 
function given in Eq. ([Tl|). We show the convergence to the asymptotic result by rewriting 
Fscat{x,y) in the form 

^scat(a;, y) = — ^ — a{x, y) , (40) 

where from Eq. (|37|) a{x, y) has the asymptotic form 

-u^{x)^f^^^'^kl'^e'^'^% + u{x) + Ud{x)- r dx'ud{x')x'V{x'). (41) 
V 4 e Jo 

The results for several values of y and the asymptotic form are shown in Fig. 3, where one 
can see that a{x, y) approaches its asymptotic form for large values of y. To better illustrate 
the convergence, in Fig. 4 we plot a(x, y) versus 1/y for x = 3. As l/y goes to zero the plot 
approaches its asymptotic value a(3, oo) = 6.863 — 4.005z given by Eq. (|l|). 

Let us now regard the asymptotic form of Fscat(a;, y) for both x and y going to infinity at 
a certain fixed angle 9 in the first quadrant. For this case only the breakup part contributes. 
Its first term F^H^ receives contributions from a saddle point and from the two end points. 
The result is 

Fscat(x, y) ^ [A{e) + B{e) + ■■■], (42) 



where 9 and p are defined by a; = pcos9,y = y 3/4 p sin 6', and 

A{9) = -he'^/^sj\ sm9e'^^''°''''^^T{kocos9) , (43) 



1 



2i 



1 



4 



B{9) = - -A{9) + A"{9) . (44) 



In the following discussion we use (x, y) and (p, 9) interchangeably. The end-point contri- 
butions are of 0{p^'^) and are beyond what is displayed in Eq. (^21). The relation between 
B{9) and A{9) is the same as found for the free propagation case considered in I. A tedious 
analytical study reveals that for this case the only contributions up to order p~^/^ are from 
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-^scat- Thus the two leading terms in Eq. (|4^) resuh solely from F^^l^ in Eq. (|20|) . The same 
is of course true for the free propagator go studied in I. 

It is clear that the first term alone in Eq. (^2]), the flux-conserving breakup behavior, is 
not a valid representation of -Fgcat at small p values. The correction term is suppressed only 
by 0[{kQp)~^], and depending upon the size of B{9) relative to A{6) the value of p may have 
to be very large before one can neglect the second- and higher-order terms in Eq. (^). To 
illustrate this property we numerically evaluate Eq. (^Tj) at pm and pm ± 10 fm for a fixed 
value of 9, which is then fit to the function 



Om H h 



(45) 



The ttin for three values of 9 along with the asymptotic value ^4(6') are given in Table I. 

From Table I one can see that Eq. (^) provides an accurate approximation to -Fscattl^;, y) 
at reasonable values of p. In addition, we note that the value of pm required for convergence 
increases as 9 increases. This feature is due to the property that large values of 9 correspond 



to small values of x, and for x small one must use Eq. (p?!). From the Taylor-series expansion 
in X of Eq. ( p7D one finds that for x greater than the range of the bound state and (x/p)^ 
small 



{k,py/^ 

For 9 = 80° the value of p must be larger than 200 fm for this approximation to be valid. 

III. THE PARTIAL-DIFFERENTIAL EQUATION APPROACH 

To test the accuracy of solving the differential form of the Faddeev equations in config- 
uration space we solve the partial differential equation 

g-'F = Q , (46) 

which has the explicit form 
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1/^2 Id Id' 



i2 



+ V{pcos9) - E 



F{p,e) = -Q{p,e), (47) 



m \dp'^ pdp p^ 86"^ 

and impose boundary conditions along a quarter circle in the first x-y quadrant at p = pmax- 
To solve the partial differential equation, we write 

F{p,e)^-^^^^^f{p,e) (48) 

and solve the resulting partial differential equation for /(p, 6). The constant /3 is a parameter 
introduced to avoid singular behavior at the origin. We solved this equation using the spline 
expansion methods described in for p < p^ax with the boundary condition that F(p, 9) 
have the form at pmax specified by Eq. (^). 

To show that this procedure can be used to obtain accurate results for the scattering 
wave function, we solved the partial differential equation numerically for various values of 
Pmax and compared the results to those obtained by numerically integrating the Green's 
function integral in Eq. (pT]). We found that values of pmax on the order of 100 fm can be 
used to obtain good wave functions. Larger values of Pmax yield more accurate solutions. In 
Figs. 5 and 6 we have plotted examples of the comparisons for fixed p values of 25 fm and 
50 fm. One can see that in both cases the agreement is excellent. 

For the three-body scattering problem one wants to obtain accurate values of the breakup 
amplitude A{9). Since this corresponds to the amplitude of the scattering wave function at 
infinity, its value cannot be obtained by evaluating the Faddeev amplitude at large values 
of Pmax- In Fig. 7 we compare the "exact" result for A{9) evaluated using Eq. (^31) and 
the integral form for T{k) given in Eq. (^) to the values extracted from the wave function 
evaluated at pmax = 200 fm. To demonstrate again that the boundary conditions for the 
partial differential equation have been treated correctly, we show the A{9) extracted from 
the numerical solution of the Faddeev equation and the evaluation of the Green's function 
integral with p = Pmax- While the wave function results are similar to the "exact" values, 
one can see that for large values of 9 they are different for the reasons discussed in the 
previous section. The A{9) obtained from the wave function at Pmax still exhibits the small 
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X behavior of the wave function for 6 near 90°. Thus, to obtain accurate results for A{6) 
one must use the integral expression. This was the procedure followed in Ref. where it 
was shown that configuration-space Faddeev calculations gave results in excellent agreement 
with the momentum-space calculations. 

IV. CONCLUSIONS 

Using a model source that mimics the real source term in the configuration-space Faddeev 
equations, the validity of the expressions for the asymptotic behavior of the wave function 
has been verified by numerically integrating the integral representation of the Green's func- 
tion integral. Using the asymptotic expressions as the boundary conditions for the partial 
differential equation in configuration space, it is possible to obtain an accurate solution of 
the scattering equation for reasonable values of Pmax- While larger values of pmax are required 
to obtain more accurate solutions, values on the order of 100 fm yield good solutions. The 
breakup amplitude corresponding to small values of x cannot be obtained from the wave 
function evaluated at large values of p; the T-matrix integral must be used to determine 
A{6) in this region. 
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Table I: Values of determined from Eq. (H) for ^ = 30°, ^ = 60°, and 6 = 80°. The 
values for = oo were determined from Eq. (^). 



Pm 


80° 


60° 


30° 


50.0 
90.0 
140.0 
190.0 
240.0 
290.0 
340.0 
oo 


Real 


Imag 


Real 


Imag 


Real 


Imag 


-0.2112 0.1240 
-0.1088 0.1718 
-0.0737 0.1651 
-0.0678 0.1609 
-0.0665 0.1586 
-0.0660 0.1578 
-0.0659 0.1574 
-0.0657 0.1560 


0.0595 0.2159 
0.0137 0.1921 
0.0083 0.1969 
0.0074 0.1980 
0.0070 0.1984 
0.0069 0.1986 
0.0068 0.1988 
0.0066 0.1990 


0.0173 0.0657 
0.0152 0.0642 
0.0155 0.0642 
0.0156 0.0642 
0.0157 0.0642 
0.0157 0.0642 
0.0157 0.0642 
0.0157 0.0642 
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FIGURE CAPTIONS 



Figure 1. The behavior of fd°^^{y) normahzed to its asymptotic value for y ^ oo. 

Figure 2. Comparison of the absolute value of fdijj) to its asymptotic form. 

Figure 3. The x dependence of a(x, y) defined in Eq. ( |iO| ) for fixed values of y and its 
asymptotic form. 

Figure 4. The real and imaginary parts of a(x, y) for x = 3 plotted versus 1/y. The triangles 
are the calculated values and the solid line is a fit to a polynomial in 1/y. 

Figure 5. Comparison of the f{p,0) evaluated using the Green's function integral and the 
/(p, 6) obtained from solving the partial differential equation for p = 25 fm. 

Figure 6. Same as for Fig. 5 for p = 50 fm. 

Figure 7. Comparison of the breakup amplitude, A{6), evaluated using the integral form 
for T{k) with the A{6) extracted from the wave function at Pmax = 200 fm. Wave function 
results for both the numerical solution of the Faddeev equation and the Green's function 
integral are shown. 
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